Kinematic power corrections in off-forward hard reactions 
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We develop a general approach to the calculation of kinematic corrections ~ t/Q^ /Q^ in hard 
processes which involve momentum transfer from the initial to the final hadron state. As the prin- 
cipal result, the complete expression is derived for the time-ordered product of two electromagnetic 
currents that includes all kinematic corrections to twist-four accuracy. The results are immediately 
applicable e.g. to the studies of deeply-virtual Compton scattering. 
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It is generally accepted that hard exclusive scattering 
processes with nonzero momentum transfer to the tar- 
get can provide one with a three-dimensional picture of 
the proton in longitudinal and transverse plane, encoded 
in generalized parton distributions (GPDs) [JQl- One of 
the principal reactions in this context is Compton scatter- 
ing with one real and one highly-virtual photon (DVCS) 
which has received a lot of attention. The QCD descrip- 
tion of DVCS is based on the operator product expan- 
sion (OPE) of the time-ordered product of two electro- 
magnetic currents where the GPDs appear as operator 
matrix elements and the coefficient functions can be cal- 
culated perturbatively. In order to unravel the transverse 
proton structure one is interested in particular in the de- 
pendence of the amplitude on the momentum transfer to 
the target t — {P' — P)^ in a reasonably broad range. 
Since, on the other hand, the available photon virtuali- 
ties are not very large, corrections of the type ~ t/Q^ 
which are formally twist-four effects, can be significant 
and should be taken into account. 

Such corrections can be called "kinematic" since they 
only involve ratios of kinematic variables and are seem- 
ingly disconnected from nonperturbative effects (e.g. one 
may consider a theoretical limit Aq^j-, ^ t ^ Q^)- 
Yet the separation of kinematic corrections ^ t/Q^ from 
generic twist-four corrections ~ Aq^q/Q^ proves to be 
surprisingly difficult. The problem is well known and 
important for phenomenology, as acknowledged by many 
authors Wld\. 

The difficulty is due to the fact that, unlike target mass 
corrections in inclusive reactions [ll| which are deter- 
mined solely by the contributions of leading twist oper- 
ators, the ~ t/Q^ corrections (and for spin-1/2 targets 
also ~ vr? IQ"^ corrections) also arise from higher- twist- 
four operators that can be reduced to total derivatives of 
the twist-two ones. Indeed, let O^^...^^ be a multiplica- 
tively renormalizable (conformal) local twist-two opera- 
tor, symmetrized and traceless over all indices. Then the 
operators 



(1) 



are, on the one hand, twist-four, and on the other hand 
their matrix elements are obviously given by the re- 



duced leading twist matrix elements, times the momen- 
tum transfer squared (up to, possibly, target mass correc- 
tions). Thus, contributions of the both operators must 
be taken into account. 

The problem arises because O2 has very peculiar prop- 
erties: divergence of a conformal operator vanishes in 
a free theory (the Ferrara-Grillo-Parisi-Gatto theorem 
[1^). As a consequence, using QCD equations of motion 
(EOM) O2 can be expressed as a sum of contributions of 
quark-antiquark-gluon operators. The simplest example 
of such relation is I13l-[l5ll 



(2) 



where O^^ = (i/2)[g7p Bv q + {fi ^ i^)] is the quark 
part of the energy-momentum tensor. The operator on 
the r.h.s. of Eq. ^ involves the gluon field strength 
and, not knowing this identity, it would be tempting to 
assume that hadronic matrix elements of this operator 
are of the order of Aq(-,j-|, which is wrong. More com- 
plicated examples involving leading-twist operators with 
two derivatives can be found in [la . [l7j . 

The general structure of such relations is, schemati- 
cally 



k 



(3) 



where Gjqk are twist-four quark-antiquark-gluon (and 
more complicated) operators and a^^'^ the numerical co- 
efficients. The subscript N stands for the number of 
derivatives in Ojq and the summation goes over all con- 
tributing operators, with and without total derivatives 
(so that in reality k is some multi- index). The same op- 
erators, Gmu, also appear in the OPE for the current 
product at the twist-four level: 



T{3{x)mY='^Y.^N,k{x)G 



Nk 



(4) 



N,k 



A separation of "kinematic" and "dynamical" contribu- 
tions implies rewriting this expansion in such a way that 
the contribution of the particular combination appearing 



2 



in (j3]) is separated from the remaining twist-four contri- 
butions. The "kinematic" approximation would corre- 
spond to taking into account this term only, and neglect- 
ing contributions of "genuine" quark-gluon operators. 

Our starting observation is that the separation of kine- 
matic and dynamical effects is only meaningful if they 
have autonomous scale dependence. Different twist-four 
operators of the same dimension mix with each other 
and satisfy a matrix renormalization group (RG) equa- 
tion which can be solved, at least in principle. Let GN,k 
be the set of multiplicatively renormalizable twist-four 
operators so that 

fc' 

The relation ([3]) tells us that one of the solutions of the 
RG equation is known without the calculation. Namely, 
there exists a twist-four operator with the anomalous di- 
mension equal to the anomalous dimension of the leading 
twist operator, and Eq. ([3]) presents the corresponding 
eigenvector. (For simplicity we ignore the contributions 
of O^On operators in this discussion; they do not pose a 
"problem" and can simply be taken into account.) 

Assume this special solution corresponds to k — 0, so 
GN,k=o = {dO)N and V'fc=o fc' ~ "^fe'- Inverting the matrix 
of coefficients i^^^^, we can write the expansion of an arbi- 
trary twist-four operator in terms of the multiplicatively 
renormalizable ones 

GN.k = 0S(9O)iv + E Q^-^' ■ (6) 

fc'#0 

Inserting this expansion in Eq. ^ one obtains 

T{j{x)mY^'' = E Civ.fc(2:)^i^o) {dO)N + . . . , (7) 

AT.fc 

where ellipses stand for the contributions of "genuine" 
twist-four operators (with different anomalous dimen- 
sions). The problem with this (formal) solution is that 

finding the coefficients cf)^^^ in general requires the knowl- 
edge of the full matrix V'i^' i alias explicit solution of the 
twist-four RG equations, which is not available. 

Twist-four operators in QCD can be divided in two 
classes: quasipartonic [l^, that only involve "plus" com- 
ponents of the fields, and non-quasipartonic which also 
include "minus" light-cone projections. Our next obser- 
vation is that quasipartonic operators are irrelevant for 
the present discussion since they have autonomous evo- 
lution (to the one-loop accuracy). Hence terms in {dO)^ 
do not appear in the re-expansion of quasipartonic opera- 
tors in multiplicatively renormalizable operators, Eq. ([S]): 
the corresponding coefficients 4>^^q vanish. As the result, 
the kinematic power correction ~ (dO)^ is entirely due 
to contributions of non-quasipartonic operators. 

Renormalization of twist-four non-quasipartonic oper- 
ators was studied systematically in |19l . |20| . The main 



result is that the RG equations can be written in terms of 
several S'L(2)-invariant kernels. Using this technique, we 
are able to prove that the anomalous dimension matrix 
for non-quasipartonic operators is hermitian with respect 
to a certain scalar product, which implies that different 
eigenvectors are mutually orthogonal, i.e. 

Z^Mfc Vl,k Vm,k ~ , (»j 

k 

where Mfe^^ is the corresponding (nontrivial) measure. 
Using this orthogonality relation and the expression ([3]) 
for the relevant eigenvector, one obtains, for the non- 
quasipartonic operators 

« = 4"^ll«^"^r^ (9) 

where jja^^^lp = J2k l^k^\'^k^^)'^ ■ Inserting this expres- 
sion in ([7]) one obtains the desired separation of kinematic 
effects. 

The actual derivation proves to be rather involved. It 
is done using the two-component spinor formalism in in- 
termediate steps and requires some specific techniques of 
the SL{2) representation theory. The purpose of this let- 
ter is to present our main result; the technical details will 
be given elsewhere. 

We define nonlocal (light-ray) vector Oy and axial- 
vector Oa operators of the leading-twist-two as the gen- 
erating functions for local twist-two operators 

0[zix, Z2x) ^[q{zix)3f: [lb] Q q{z2x)\^ ^ . (10) 

Here is a four- vector which is not necessarily light-like, 
zi and Z2 are real numbers and Q is the (diagonal) matrix 
of quark electromagnetic charges squared. The Wilson 
line between the quark fields is implied. The leading- 
twist projector [. . .];.t. stands for the subtraction of traces 
of the local operators so that formally 

[q{zix):l:Qq{z2x)\^^^^ = 

N 

+Z2i:)p„]Q(?(0) - tracesj. (11) 

The leading-twist light-ray operators satisfy the Laplace 
equation d'^0{zix, Z2x) — . Explicit form of the projec- 
tor [. . .];.(. is irrelevant for the further discussion (some 
useful representations can be found e.g. in (gI [2^). 

Alternatively, one can expand a nonlocal operator in 
the contributions of multiplicatively renormalizable (in 
one loop) conformal operators 

0{ziX, Z2X) ^'S^ Kpf Z^2 / du{uu)^'^^[ON{Zi2x)]^ ^ , 

(12) 

where xn = 2{2N + 3)/(Ar -1-1)!. Here and below we use 
the following shorthand notation: 

U = 1 — U , Zi2 — Zi — Z2 , Z^2 — + UZ2 . (13) 
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The conformal operator Oat is defined as 



X 0{zix + y, z-ix + 2;) 



(14) 



where C^^{x) is the Gegenbauer polynomial. 

We are able to find the contributions related to the 
leading-twist operator ()10p in the time-ordered product 
of two electromagnetic currents 



^T{f;"{x)jr{0)} (15) 
to the twist-four accuracy. The resuh has the form 



TT a; 1^ L 

{x^,d^ + x^d^,)X + {x^d^ ~ x„dfj,)Y | , (16) 



where 9^ = d/dx^", S^ai^p = gtiagup + gi^ag^ip - g^Lugaf} 

and the totally antisymmetric tensor is defined such that 
£0123 = 1- The expansion of and starts from twist 
two, = V*f 2 + V^=3 + V^=4 ^ . . ., A/j 



-I- 



while 



. . . , r^p — + A^ 

and Y = Y*='^ + . . . 



twist-four. 

It turns out that vector operators always appear to be 
antisymmetrized and axial- vector symmetrized over the 
quark and antiquark positions, respectively, so we define 
the corresponding combinations: 

0'f\zi,Z2) = [q{zix)^Qq{z2x)]^^^ - {zi ■(^ Z2) , (17) 
0^j^\zi,Z2) = [q{zix):^j5Qq{z2x)]^^+{zi^Z2). 

The leading-twist expressions are well known and can be 
written as (cf. 22]) 



t=2 



1 



da 



^.^j^ duO\,\u,0), 



(18) 



Note that separation of the leading-twist contributions 
[. . .]i,t. from the nonlocal operators by itself produces a 
series of kinematic power corrections to the amplitudes, 
which are analogous to Nachtmann target mass correc- 
tions to deep- inelastic scattering [llj. Such corrections 
are discussed in detail in 

For the twist-three functions we obtain 



Tt = 3 



\-aud^'x^d'']0 



[u)]] . (19) 



Here Pj/ is the momentum operator 

[zP.,q(y)] = {p'mM\p) = ip' -p)Ap'\0\p) 



and we used a notation 



0\^\z)^\j\wO(^\z,w). 



(20) 



One can easily check that x^V^ 
similarly x^'h}=^ = a^A^f^ ^ q. 



and 

Note that the terms 
in Init in Eqs. (jl9p are themselves twist-four and can be 
dropped if the calculation is done to twist-three accuracy. 
The resulting simplified expression is in agreement with 
Refs. PiQ- These terms must be included, however, in 
order to ensure the separation of twist-three and twist- 
four contributions. 

The twist- four contributions V^^*, A^^^, X*^^ and 
Yt=4 present our main result. In this case we did not 
find a simple nonlocal representation and write the an- 
swer in terms of integrals over the position of the local 
operators, cf. Eq. ([T2|) . This form is equally well known 
and usually referred to as conformal OPE [231. It proves 
to be the most convenient for implementing the scale de- 
pendence in leading- twist GPDs . We obtain 



V!: 



Y* 



2 ^ "" (N + 2Y Jo 



du{uuf+Ux^[O^Niux)]u. +i-N{N + S) I dvv''-'x'd^[O^Niuvx)]i.t. S> , 



.N-1 2' 



4e 



xn- 



N + 1 



1 nl 



A ^ ""('/V-^2^2 / {u-u) \ dvv [0,,j{uvx)]i,t. , 

AT.odd \ + ) Jo Jo 



AT.odd 



1 fl 

Ni,.2 , -2^ / j„,„,7V-l \rnV I 



''^TW^ I du{uu)^iu' + u') I dvv^-'[0^^{uvx)]i.t. 
(N + 2y Jo Jo 



(21) 
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Here On is defined as the divergence of the leading-twist 
conformal operator, cf. O2 in Eq. ([T|): 



ON{y) 



1 



d 



(22) 



One of the two integrals over the position of Oat can eas- 
ily be taken, resulting in slightly more lengthy expres- 
sions. 

Note that the operator Oi in Eq. ([T]) , which reads 
[iP ^[iP^ ,On] in our present notation, does not con- 
tribute to the answer for our special choice of the cor- 
relation function T{j^[x)jv{0)}. The T-product with 
generic positions of the currents, T{j^{zix)jv{z2x)}, in- 
cludes both operators. The corresponding result is much 
more cumbersome and will be given elsewhere. 

For comparison we rewrite the leading-twist contribu- 
tion in the same form: 



V, 



^ N- 

N,odd 



duu^u^+^Oliux), (23) 



where C']^(mx) is the conformal operator (|14p at the 
space-time position ux. 

Conservation of the electromagnetic current implies 
that 

d^T^.ix) - , d-'T^^ix) - ^[P^T^,(a;)] . (24) 

We have checked that these identities are satisfied up to 
twist- 5 terms. 

For completeness we give the expression for the oper- 
ator [iP^, d^O{zi, Z2)] entering the twist-three functions 



Yt-3 



vt-3 



in terms of O 



N- 



1 

[zP,„a''0(zi,Z2)] - -5+/urfu[iP^[iP^,0(uZl,MZ2)]] 

+ ^KN{N+lfz^2 I dvv^ [ du {uu)^+^dN{vz'l2x), 
Jo Jo 

(25) 

where 5"+ = zfd^-^ + + 2zi + 2z2. In phenomenolog- 
ical applications it can be advantageous to use relations 
of this kind to rewrite all contributions of Oat in terms 

of [iP^,9''0(zi,Z2)]. 

To summarize, we have given a complete expression 
for the time-ordered product of the two electromagnetic 
currents that resums all kinematic corrections to the 
twist-four accuracy. The results have immediate applica- 
tions to the studies of deeply- virtual Compton scattering 
and 7* — >■ (tt, ?7, . . .) -f- 7 transition form factors. The 
twist-four terms calculated in this work give rise both 
to a t/Q^ correction and the target mass correction 
^ m? /Q^ for DVCS, whereas for the transition form fac- 
tors these two effects are indistinguishable as there is only 
one mass scale. We remark that the distinction between 
the kinematic corrections due to contributions of leading- 
twist [4-10] and higher-twist operators considered in our 
work is not invariant under translations along the line 
connecting the currents and has no physical meaning. 
Such corrections must always be summed up. Concrete 
applications go beyond the tasks of this letter. 
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